clear
 use "$P_Data/GS/gas_station.dta"
 by Brand, sort: gener n_stations = _N
by Brand, sort: keep if _n==1
 sort n_stations
 keep if _n>(_N-5)
  keep Brand n_stations
  save "$P_Data_Processed/top5_brands.dta", replace
